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ABSTRACT 

The  general  characteristics  and  the  bandwidth  response 
of  the  Kaijo  Denki  PAT  313-1  sonic  anemometer / thermomet er 
were  examined.   Measurements  were  compared  to  those  of  a 
hot  wire  anemometer  system  and  a  one  mil  thermocouple  sys- 
tem.  Turbulence  data  were  obtained  by  all  three  systems 


sim 


ultaneously  in  two  experiments  conducted  at  a  beach 


site.   The  sonic  anemometer /thermomet er  system  proved  to 
be  a  reliable,  rugged,  easily  installed  system  with  a  band^ 
width  response  of  from  0.0  Hz  to  approximately  12.0  Hz. 


TABLE    OF    CONTENTS 

I.  INTRODUCTION 11 

A.  PROBLEM 11 

B.  OBJECTIVES  11 

II.  GENERAL  DESCRIPTION  AND  OPERATING  PROCEDURES 
OF  THE  SONIC  ANEMOMETER/THERMOMETER  MODEL 

PAT  313-1 12 

A.  GENERAL 12 

B.  COMPONENT  PARTS  OF  THE  SONIC 
ANEMOMETER/THERMOMETER  12 

1.  Power  Supply/Indicator  Unit  (PR-1)  12 

2.  Probe  (TR-31)  13 

3.  Junction  Box  (OA-33H)  18 

C.  ADVANTAGES/DISADVANTAGES  OF  THE  SONIC 
ANEMOMETER/THERMOMETER  18 

D.  MEASUREMENT  PROCESS  24 

1.  Wind  Measurement  24 

2.  Temperature  Measurement  26 

3.  Discussion  of  Possible  Errors  in 
Wind  Velocity  and  Temperature 
Measurements  27 

4.  Timing  Sequence  30 

III.  DATA  ANALYSIS  34 

A.  GENERAL  34 

B.  KOLMOGOROV'S  HYPOTHESES  35 

1.  General  35 

2.  Kolmogorov's  First  Hypothesis: 
Equilibrium  Range  36 


3.   Kolmogorov's  Second  Hypothesis: 

Inertial  Subrange  36 

C.  ANALOG  TO  DIGITAL  PROCEDURES  37 

1.  General  37 

2.  Considerations  Prior  to  Digitization  39 

D.  CONVERSION  OF  SEVEN  TRACK  TO 

NINE  TRACK 41 

E.  FOURIER  COEFFICIENT  COMPUTATIONS  42 

F.  UBCSCOR  PROGRAM 43 

IV.  FIELD  EXPERIMENTS  45 

A.  PRELIMINARIES  45 

B.  OBJECTIVE  OF  EXPERIMENTS  45 

C.  EXPERIMENTS 46 

1.  Installation  of  Instruments  46 

2.  Experiment  1-14  Dec  1972  46 

3.  Experiment  II-  27  Feb  1973 51 

V.  RESULTS  AND  CONCLUSIONS  54 

A.  GENERAL  PERFORMANCE  OF  SONIC 
ANEMOMETER/THERMOMETER 54 

B.  SPECTRA  PLOTS  54 

1.  General 54 

2.  Sonic  Anemometer /Thermometer  Data 
Sampled  at  390.5  Hz  -  14  Dec  1972 
(Experiment  I)  56 

3.  Sonic  Anemometer /Thermometer  Data 
Sampled  at  62.5  Hz  -  14  Dec  1972 
(Experiment  I)  56 

4.  Hot  Wire  Anemometer  Data  Sampled 

at  200.0  Hz  -  10  Jan  1973  56 


5.  Sonic  Anemometer /Thermometer  Data 
Sampled  at  62.5  Hz  -  28  Feb  1973 
(Experiment  II)  79 

6.  Hot  Wire  Anemometer  Data  Sampled 
at  200.0  Hz  -  28  Feb  1973 

(Experiment  II)  79 

7.  Thermocouple  Data  Sampled  at  200.0  Hz 

-  28  Feb  1973  (Experiment  II)  79 

C.   CONCLUSIONS  79 

APPENDIX  A  -  Definition  of  Terms  Used  in 

Analog-Digital  Data  Processing  94 

BIBLIOGRAPHY  96 

INITIAL  DISTRIBUTION  LIST  97 

FORM  DD  1473  98 


LIST  OF  FIGURES 

1.  Wind  and  Sonic  Axes  14 

2.  Power  Supply /Indicator  Unit  (PR-1)  15 

3.  Wind  Unit  (PD-11)  16 

4.  Temperature  Unit  (PS-11)  17 

5.  Wind  Direction  Limits  of  Probe  Type 

TR-31.  Kaijo  Denki  (1969)  19 

6.  Probe  Head.   Kaijo  Denki  (1969)  19 

7.  Probe  Type,  TR-31.   Kaijo  Denki  (1969)  20 

8.  Junction  Box  (OA-33H) .  Kaijo  Denki  (1969)  21 

9.  Wiring  Diagram  of  PAT-313-1  System. 

Kaijo  Denki  (1969)  21 

10.   Schematic  of  Typical  Measuring  Process 

and  Nomenclature  25 

11A .  Temperature  Compensating  Constant. 

Kaijo  Denki  (1969)  29 

11B.  Humidity  Compensation  Graph. 

Kaijo  Denki  (1969)  29 

12.  Error  by  Wind  Velocity. 

Kaijo  Denki  (1969)  30 

13.  Timing  Sequence  Schematic.   Kaijo  Denki  (1969)  32 

14.  Temperature  and  Time  Relationship. 

Kaijo  Denki  (1969)  33 

15.  Schematic  Drawing  of  Spectral  Energy 
Transfer  in  Free  Convection.  Lumley  and 

Panofsky  (1964)  38 

16.  Arrangement  of  Probes  for 

Experiments  I  and  II  47 

17.  Close-up  of  Probe  Arrangement  48 


18. 
19. 
20. 

21. 
22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 
30. 
31. 


Probe  Stand  Orientation  49 

Inside  Installation  of  Equipment  50 

Strip  Chart  Recording  of  the  Sonic  A 

and  W  Channels  Indicating  Spikes  52 

Sonic  "A"  Channel  Data  Observed  14  December 

1972  1st  Period  (Experiment  I) 

High  Sample  Rate  57 

Sonic  "B"  Channel  Data  Observed  14  December 

1972  1st  Period  (Experiment  I) 

High  Sample  Rate  58 

Sonic  "W"  Channel  Data  Observed  14  December 

1972  1st  Period  (Experiment  I) 

High  Sample  Rate 59 

Sonic  "T"  Channel  Data  Observed  14  December 

1972  1st  Period  (Experiment  I) 

High  Sample  Rate 60 

Sonic  "A"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

High  Sample  Rate  61 

Sonic  "B"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

High  Sample  Rate  62 

Sonic  "W"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

High  Sample  Rate 63 

Sonic  "T"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

High  Sample  Rate  64 

Sonic  "A"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

High  Sample  Rate  65 

Sonic  "B"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

High  Sample  Rate  66 

Sonic  "W"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

High  Sample  Rate  67 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


Sonic  "T"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

High  Sample  Rate  68 

Sonic  "A"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

Low  Sample  Rate  69 

Sonic  "B"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

Low  Sample  Rate  70 

Sonic  "W"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

Low  Sample  Rate  71 

Sonic  "T"  Channel  Data  Observed  14  December 

1972  2nd  Period  (Experiment  I) 

Low  Sample  Rate 72 

Sonic  "A"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

Low  Sample  Rate  73 

Sonic  "B"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

Low  Sample  Rate  74 

Sonic  "W"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

Low  Sample  Rate  75 

Sonic  "T"  Channel  Data  Observed  14  December 

1972  3rd  Period  (Experiment  I) 

Low  Sample  Rate  76 

Hot  Wire  Data  Observed  10  January  1973 

1st  Period  77 

Hot  Wire  Data  Observed  10  January  1973 

2nd  Period  78 

Sonic  "A"  Channel  Data  Observed  28  February 

1973  1st  Period  (Experiment  II) 80 

Sonic  "W"  Channel  Data  Observed  28  February 

1973  1st  Period  (Experiment  II)  81 

Sonic  "T"  Channel  Data  Observed  28  February 

1973  1st  Period  (Experiment  II)  82 


46. 


47. 


48. 


49 


50 


51. 


52 


53 


54 


55 


Sonic  "A"  Channel  Data  Observed  28  February 

1973  2nd  Period  (Experiment  II)  83 

Sonic  "W"  Channel  Data  Observed  28  February 

1973  2nd  Period  (Experiment  II)  84 

Sonic  "T"  Channel  Data  Observed  28  February 

1973  2nd  Period  (Experiment  II)  85 

Sonic  "A"  Channel  Data  Observed  28  February 

1973  3rd  Period  (Experiment  II)  86 

Sonic  "W"  Channel  Data  Observed  28  February 

1973  3rd  Period  (Experiment  II)  87 

Sonic  "T"  Channel  Data  Observed  28  February 

1973  3rd  Period  (Experiment  II)  88 

Uot  Wire  Data  Observed  28  February  1973 

2nd  Period  (Experiment  II)  89 

Hot  Wire  Data  Observed  28  February  1973 

3rd  Period  (Experiment  II)  90 

Thermocouple  Data  Observed  28  February 

1973  2nd  Period  (Experiment  II)  91 

Thermocouple  Data  Observed  28  February 

1973  3rd  Period  (Experiment  II)  92 


ACKNOWLEDGEMENT 

The  author  acknowledges  the  generous  assistance  and 
advice  given  by  Dr.  Kenneth  L.  Davidson  and  Dr.  Noel  E.J. 
Boston  on  all  phases  of  this  research.   A  special  note  of 
thanks  to  Dr.  Thomas  M.  Houlihan  for  his  expertise  and 
assistance  with  the  hot  wire  anemometer  and  one  mil  ther- 
mocouple systems  during  the  field  experiments,  Mr.  Robert 
Limes  for  his  assistance  during  the  digitization  of  data, 
Miss  Sharon  Raney  for  advice  in  IBM  tape  procedures  and 
FT1  Jim  Nolan  USN  for  his  expert  maintenance  and  upkeep 
of  the  sonic  anemome t er / thermome t er . 

This  research  project  was  supported  by  the  Naval  Ord- 
nance Systems  Command  through  project  order-3-0126  dated 
1A  September  1972  with  the  Naval  Postgraduate  School. 


10 


I.    INTRODUCTION 

A.  PROBLEM 

Turbulent  processes  in  the  atmosphere  are  of  interest 
to  meteorologists  since  it  is  by  these  processes  that  effi- 
cient vertical  and  horizontal  transfers  of  energy,  heat 
and  momentum  occur.   These  transfer  processes  have  been 
experimentally  observed  to  occur  within  the  frequency  band 
from  0.01  Hz  to  10.0  Hz.   An  important  method  for  describ- 
ing the  role  of  turbulence  is  through  use  of  a  system  of 
instruments  capable  of  measuring  over  this  band.   The  ana- 
log data  obtained  is  then  read  into  computer  systems  to  be 
digitized  into  discrete  samples  and  then  to  be  analyzed  by 
spectral  or  other  techniques. 

B.  OBJECTIVES 

The  primary  objective   of  this  thesis  was  to  utilize  a 
sonic  anemome t er / thermome ter  to  measure  turbulence  and 
evaluate  the  band  over  which  it  is  capable  of  measuring. 
A  secondary  objective  was  to  utilize  a  Hybrid  computer  sys- 
tem to  digitize  the  analog  data  and  then  to  use  the  IBM- 
360/67  computer  system  for  spectral  analysis  of  the  digi- 
tized data.   A  third  objective  was  to  compare  a  sonic 
anemome te r/ the rmome t er  with  a  hot  wire  anemometer  system 
and  a  one  mil  thermocouple  system. 
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I I .   GENERAL  DESCRIPTION  AND  OPERATING  PROCEDURES 

OF  THE 
SONIC  ANEMOMETER/THERMOMETER  MODEL  PAT-313-1 


A.  GENERAL 

The  sonic  anemometer / thermomet er  "  measures  the  three 
velocity  components  of  the  wind  by  detecting  accurately 
the  transit  time  d  if f erence  (transit  time  sum  for  tempera- 
ture) between  two  acoustic  pulses  propagated  in  opposite 
directions  across  a  fixed  path.   The  three  velocity  com- 
ponents of  the  wind  are  the  longitudinal,  lateral  and  ver- 
tical components  denoted  by  u,  v,  and  w  respectively  alonj 
axes  X,  Y  and  Z  (Figure  1) .   The  w  component  is  obtained 
directly  from  the  vertical  axis,  but  u  and  v  have  to  be 

computed  from  the  velocity  components  U,  and  V_  measured 

A        B 

along  the  horizontal  sonic  paths  A  and  B  which  are  120 
degrees  apart . 

B.  COMPONENT  PARTS  OF  THE  SONIC  ANEMOMETER/THERMOMETER 
1  .   Power  Supply/Indicator  Unit  (PR-1) 

The  power  supply/ indicator  unit  and  its  component 
breakdown  is  pictured  in  Figures  2,  3,  4.   The  instrument 
operates  on  100/115  volts  with  +  10%  variation  at  50  to 
60  Hz.   Total  power  required  is  100  watts.   All  switches 
and  controls  are  located  in  this  unit.   Calibration  and 
adjustment  procedures  (discussed  later)  are  performed  at 


Model  PAT-313-1  manufactured  by  the  Kaijo  Denki  Company 
Limited  of  Tokyo,  Japan. 
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this  unit.   Visual  meter  output  (analog  form)  is  available 
for  wind  measurement  on  three  range  scales  of  1)  0  to  +  1 
m/sec,  2)  0  to  +  3  m/sec  and  3)  0  to  +  10  m/sec       (Fig- 
ure 3).   Temperature  meter  output  is  +  IOC  from  mean  air 
temperature  which  is  dialed  in  on  8  dial  (Figure  4).   Two 
BNC  output  jacks  are  available  for  tape  recording  or  sig- 
nal monitoring  purposes  on  both  wind  and  temperature  units. 
2  .   Probe  (TR-31) 

The  probe  has  to  be  employed  within  the  wind  direc- 
tion limits  indicated  in  Figure  5.   The  wind  velocity  would 
be  less  than  actual  if  the  probe  were  positioned  so  that 
the  wind  was  parallel  to  the  path  travelled  by  the  sound 
pulses  or  if  the  supporting  structure  of  the  probe  inter- 
rupted the  wind  flow. 

The  probe  heads  are  ceramic  type  transducers  from 
which  sound  waves  are  trasmitted  at  a  frequency  of  lOOKKs. 
The  transducers  have  high  directivity  and  damping  to  mini- 
mize interference  with  other  channels  and  are  small  physi- 
cally in  order  to  avoid  blocking  the  air  flow  (Figure  6). 
The  transducers  screw  into  special  holders  designed  to 
minimize  acoustic  coupling  through  the  array.   All  probe 
heads  are  interchangeable  with  only  minor  adjustments  re- 
quired in  the  wind/ temper at ure   units  affected  by  the 
change . 

The  probe  heads,  holders  and  entire  array  structure 
have  been  designed  to  reduce  interference  to  the  wind  flow. 
All  wiring  and  connections  have  been  inserted  inside  the 
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Figure    1.       Wind    and    Sonic    Axes. 
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Figure  2.   Power  Supp ly /Indicator  Unit  (PR-1) . 
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Figure    3.       Wind    Unit     (PD-11)  . 
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probe  tubing  for  protection  from  the  elements  as  well  as 
for  reducing  air  flow  interference  (Figure  7). 
3 .   Junction  Box  (OA-33H) 

The  junction  box  contains  three  sets  of  pulse  os- 
cillators (which  generate  the  lOOKHz  signal  sent  to  the 
transducers)  and  six  sets  of  head  amplifiers  (which  amplify 
by  1000  the  lOOKHz  signal  sent  to  the  transducers)  .   The 
junction  box  is  a  watertight  unit  with  connectors  at  both 
ends  (Figure  8).   One  hundred  meters  of  cable  link  the 
junction  box  and  supply / indicator  unit  and  three  meters  of 
cable  link  the  junction  box  and  probe  (Figure  9).   The 
junction  box  mus  t  be  connected  to  the  probe  and  power  sup- 
ply prior  to  turning  the  sonic  anemometer  on  or  the  output 
transformer  will  be  shorted. 

C.   ADVANTAGES/DISADVANTAGES  OF  THE  SONIC  ANEMOMETER/ 
THERMOMETER 

The  sonic  anemome ter/ thermometer  has  many  advantages 
over  previously  developed  wind  and  temperature  measuring 
systems.   Recent  results  by  McBean  (1972)  revealed  the 
advantages  the  sonic  anemometer  has  over  a  Gill  propeller 
anemometer  especially  with  respect  to  sensor  response  and 
angle  of  -attack  (cosine  response)  in  eddy  correlation 
measurements.   A  report  by  Kaimal  e_t_  al_  (1971),  in  which 
they  compared  the  pulse  and  continuous-wave  type  sonic 
anemometers,  lists  many  advantages  of  the  sonic  anemo- 
meter system  over  other  velocity  measuring  systems. 
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Figure  5.   Wind  Direction  Limits  of  Probe  Type  TR-31. 
Kaijo  Denki  (1969) 
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Figure  6.   Probe  Head.   Kaijo  Denki  (1969) 
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Fieure  7.   Probe  Type  TR-31.   Kaijo  Denki  (1969) 
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Figure  8.   Junction  Box  (OA-33H). 
Kaijo  Denki  (1969) 
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Figure    9.       Wiring    Diagram    of    PAT -313-1    System. 
Kaijo    Denki     (1969) 
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The  PAT-313-1  system  has  no  mechanical  parts  which 
eliminates  lag  due  to  mechanical  inertia  and  thus  allows 
the  sonic  to  measure  very  light  and  sudden  changes  of  the 
wind.   The  sonic  also  has  high  sensitivity  and  a  broad 
range  in  all  three  wind  components  and  in  the  temperature 
component . 

The  resolution  of  the  sonic  anemometer  is  +  2  cm/sec 
in  a  span  of  0-10  m/sec  and  the  resolution  of  the  sonic 
thermometer  is  +  0.1C  in  a  span  of  +  IOC  centered  around 
a  mean  temperature  range  of  -IOC  to  AOC.   The  response  is 
linear  and  so  is  its  sensitivity.    (Therefore,  the  res- 
ponse is  the  same  at  5  cm/sec  as  it  is  at  10  cm/sec.) 

The  sonic  anemometer  has  a  time  response  in  the  milli- 
second range  as  it  detects  changes  in  pulse  time  between 
two  sets  of  transducers.   It  takes  measurements  440  times 
/sec  for  the  wind  components  and  220  times/sec  for  the 
temperature  component  which  allows  it  to  detect  changes  of 
approximately  100  Hz  at  an  accuracy  of  three  percent 
(Kiajo  Denki,  1969)  . 

The  sonic  anemometer  is  not  affected  by  exposure  to 
the  direct  rays  of  the  sun  or  to  the  rain  as  would  a  thermo- 
meter or  a  hot  wire  anemometer.   It  is,  however,  affected 
by  immersion  in  water.   A  high  moisture,  high  heat  environ- 
ment  affects  measurement  procedures.   These  have  to  be 
corrected  in  order  to  get  true  readings  because  the  measur- 
ing circuits  are  set  for  an  atmosphere  of  zero  percent 
relative  humidity  at  a  temperature  of  20C.   The  probe, 
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junction  box  and  connecting  cables  are  waterproofed  allow- 
ing for  all  weather  operation  of  this  system.   There  are 
no  easily  broken,  sensitive  elements  (such  as  a  hot  wire) 
exposed  to  the  weather.   It  is  a  rugged  system  that  can 
accept  some  rough  handling.   Cable  connection  points,  es- 
pecially where  they  enter  the  probe,  are  sensitive  areas. 
Too  much  twisting  of  the  cables  can  result  in  broken  wires 
in  the  connectors.   Care  must  be  exercised  so  as  not  to 
bend  the  probe  thus  changing  the  length  of  transmission  of 
the  acoustic  path  and  causing  misalignment  of  the  probe 

heads.   Both  problems  would  create  serious  errors  in  velo- 

2 
city  and  temperature  measurement.    Each  of  the  probe  heads 

has  a  rubber  cap  over  it  to  protect  the  ceramic  transducer 

element  but  abrasive  handling  will  rip  it. 

The  probe  has  a  small  bubble  level  in  the  top  plate 

which  is  used  to  insure  that  it  is  properly  leveled  when 

installed.   Errors  are  introduced  if  the  probe  is  not  hori- 

3 

zontal.    The  probe  should  be  mounted  in  such  a  way  that 

it  is  free  of  vibrations.  It  measures  relative  speed  of 
the  environment  with  respect  to  the  probe  so  any  movement 
of  the  probe  will  introduce  errors  into  the  signals.  Any 
object  placed  within  two  meters  of  the  probe  will  reflect 
erroneous  signals  into  the  transducers.  These  are  recog- 
nizable as  noise  spikes. 


'Kaimal  e_t_  al_  (1971)  has  an  excellent  discussion  of  these 
errors . 

Kaimal  et  al  (1971)  discusses  this  error. 
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D.   MEASUREMENT  PROCESS 
1 .   Wind  Measurement 

As  described  previously,  the  sonic  anemometer  mea- 
sures wind  velocity  by  detecting  the  transit  time  differ- 
ence between  two  acoustic  pulses  propagated  in  opposite 
directions  across  a  fixed  path.   The  nomenclature  used  to 
describe  the  measuring  process  for  a  typical  pair  of 
transducers  is  shown  in  Figure  10. 

The  propagation  times  (t-,,  t„)  for  the  pulses  from 
TRX  to  R   and  TR   to  R_  are 


t,  = 


1    C  cos  a  +  V 


(1) 


and 


t„  - 


2    C  cos  a  -  V 


(2) 


The  travel  time  (t-,)  of  the  sound  with  the  wind  is  less 
than  the  travel  time  (t_)  with  the  sound  against  the  wind. 
The  difference  A  t  =  t  „  -  t   is  given  by 


2 IV 


At  = 


(C  cos  a)2  -  V  2 
a 


(3) 


Angle  a  is  assumed  to  be  small;  therefore,  cos  a 

2        2 
""  1.   One  further  assumes  that  C   >  >  V   and  therefore, 

2         2 
C   >  >  V   .   These  approximations  result  in 


"7T' 


(A) 


which,  when  solved  for  V   gives 


V  .  = 


C2At 


d     2£ 


(5) 
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V    Wind  Velocity 
and 
Sound  Transmitting 


TR 


C  cos  a  t 


C  cosa  t 


KEY: 


I    =  distance  between  transmitter  and  receiver 
C  =  sound  velocity 
V  =  wind  velocity 


V   =  component  of  wind  velocity  perpendicular 
to  sound  transmitting  path 

V,  =  component  of  wind  velocity  parallel  to 
sound  transmitting  path 

t..  =  time  for  sound  to  travel  from  T ..  to  R, 
(with  wind) 

t„  =  time  for  sound  to  travel  from  T~  to  R» 
(against  wind) 

TR  ,  TR   =  transmitting  transducers 

R..  ,  R~  =  receiving  transducers 

a  =  angle  between  transmitting  path  and  dis- 
tressed sound  path 

Figure  10.   Schematic  of  Typical  Measuring  Process  and 
Nomenclature . 
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The  wind  velocity  parallel  to  the  sound  transmitting  path 

(V.)  is  therefore  determined  and  is  the  quantity  indicated 

on  the  metered  output  of  the  sonic  anemometer  as  U   or  V  . 

A       B 

The  values  of  u  and  v  (the  longitudinal  and  lateral  wind 
components  of  the  X  and  Y  axes  of  Figure  1)  have  to  be  com- 
puted from  the  velocity  components  U   and  V  .   The  w  com- 

A        B 

ponent  is  taken  directly  from  the  output. 
The  equations  for  u  and  v  are: 

u  =  (cos  $  -  -  sin  3)  UA  +  (cos  3  +  -  sin  3)  VD      (6) 
/J  A  /3  B 


v  =-(  -  cos  3  +  sin  3)  UA+  (  -  cos  3  -  sin  3)V„ 
•3  A    /3  B 


(7) 


where 


3  =  tan 


VB  -  UA 


/3  (VB  +  UA) 


Overbar  denotes  time  average  over  a  suitably  long 
period  of  time   Kaimal  e_t  al_  (1971). 
2 .   Temperature  Measurement 

Temperature  is  measured  in  a  r.snner  similar  to  the 
wind  measurement  except  that  the  propagation  times  (t, ,t.) 
are  added  instead  of  subtracted.   The  w  component  trans- 
ducers provide  the  data  used  for  this  measurement.   The  de' 


fining  expression  for  temperature  is 


t  =  t1+  t2  = 


2      2  1/2 
2A(C   -  V  z)1/z 
n_ 

2      2 

C   "  Vd 


(8) 
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2  2  2  2 

Again    assuming    C       >    >    V  and    C       >    >    V 

a  n 


equation  (8)  reduces  to 


21 
C 


(9) 


The  speed  of  sound  in  air  is  expressed  as 
C  =  20.06  •  [T  +  (1  +  0.32  -)]1/2 


(10) 


where  T  is  in  degrees  absolute  (K)  ,  e  is  the  vapor  pressure 
in  atmospheres  and  p  is  the  atmospheric  pressure  in  atmos- 
pheres.  From  (9)  and  (10) 


T  = 


25,  2 
(~) 


402.68(1  +  0. 32  -) 

P 


(ID 


Furthermore,  set  T_„  =  T(l  +  0.32  -)  where  T„„  is 

bv  p  S  V 

sound  virtual  temperature  (taken  as  the  true  temperature 
in  degrees  K)  measured  by  the  sonic  thermometer.   Equation 
(11)  is  rewritten  as: 


=  (• 


2£    1/2 


(12) 


SV    v20.067t 
T    is  the  parameter  obtained  by  measuring  t  from  the  son  c 

D  V 

thermometer.   Being  a  virtual  temperature,  this  measure- 
ment is  subject  to  humidity  effects. 

3 .   Discussion  of  Possible  Errors  in  Wind  Velocity 
and  Temperature  Measurements 

Errors  are  introduced  by  humidity,  temperature  and 
atmospheric  pressure  which  require  correcting  the  readings 
from  the  sonic  anemometer/thermometer.   Detailed  discussion 
of  errors  resulting  from  misalignment  of  the  array,  line- 
averaging  ,  f  in  ite  spatial  resolution  of  the  array  and  noise 
in  the  detector  has  been  described  by  Kaimal  e  t  al  (19  71)  . 
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As  seen  in  equation  (5)  and  equation  (10),  the  wind 
velocity  measured  will  be  a  function  of  temperature,  humid- 
ity and  atmospheric  pressure.   From  equation  (1?),  however, 
the  measured  temperature  will  depend  on  humidity  and  atmos- 
pheric pressure  only.   For  hot  humid  days  the  effect  on 
velocity  cannot  be  ignored.   Humidity  significantly  effects 
the  absolute  temperature  but  is  negligible  as  far  as  fluctua^ 
tions  in  temperature  are  concerned  provided  the  fluctuations 
in  humidity  are  small. 

The  ratio  K,  between  the  wind  velocity  V    (at  the 
temperature  T    +  293K  or  20C  and  e/p  =  0)   and  the  wind 
velocity  V     (at  the  temperature  T   )  is 


V 


K  = 


TSV 


SV 


293  +  6 


SV 


293 


(13) 


20     20 
where  0    =  sound  virtual  temperature  measured  by  the  sonic 

O  V 

thermometer  in  degrees  C.   Any  wind  velocity  or  temperature 
measured  by  the  sonic  anemome ter/ thermometer  in  an  atmos- 
phere other  than  20C  and  zero  percent  relative  humidity 

must  be  corrected  by  a  factor  K  (Figure  11A,  11B). 

2      2 
The  error  introduced  by  the  assumption  of  C   >>  V 


is 


Zw  = 


(14) 


and  is  such  that  for  wind  speeds  of  less  than  10  m/sec  it 
is  zero  (Figure  12). 


Sonic  anemome t er/ t hermome ter  circuitry  is  set  for  a  base 
temperature  of  20C  and  zero  percent  relative  humidity. 
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Figure  11A.   Temperature  Compensating  Constant.  . 
Note:   The  9  dial  and  the  meter  scale  of  the  TEMP  UNIT  calibrated  in  9SV- 
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Figure    11B.       Humidity    Compensation    Graph. 
Kaijo    Denki     (1969) 
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(The   temperature 
measured  by  ultrasonic 
thermometer) 

RH  =  Relative 
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4 .   Timing  Sequence 
A  more  detailed 
analysis  by  Soumie  and 
Businger  (1958)  showed 
that  velocity  and  tempera- 
ture fluctuations,  within 
certain  limits,  were  pro- 
portional to  the  fluctua- 
tions in  the  difference 
and  the  sum  of  the  time  in- 


iv2 


7 


1 

1 

— 

1 

10   20       30       40       SO 

V  (mAec) 


Figure  12.   Error  by  Wind 
Velocity.  Kaijo  Denki  (1969) 


tervals  respectively.  A  discussion  of  the  method  of  obtain- 
ing these  time  intervals  is  of  use  in  order  to  understand 
the  setting  up  and  calibration  of  the  sonic  anemome ter/ ther- 
mometer. Only  the  vertical  transmitter-receiver  array  (that 
used  to  measure  the  vertical  velocity  component  and  tempera- 
ture) needs  to  be  considered  in  order  to  understand  the 
triggering  and  pulse  measuring  system  used  in  the  instrument 

Pulses  are  emitted  simultaneously  from  TRTT1  and  TRTTO 

Wl         w  Z 

at  time  T,  and  reach  their  respective  receivers  RT71  and  RT7„ 

Wl       Wz 

after  time  intervals  of   t  ..  and  t  ~   seconds  (Figure  13)  . 

wl       w2 


The  time  difference     At   =  t  „  -  t  ..         is  converted 

w     w2     wl 

by  a  wind  analog  card  (within  the  w  unit)  and  changed  into 

a  voltage  representing  the  vertical  component  of  the  wind 

velocity.   The  output  of  R  ,  drives  the  delay  card  Dn  which 

wl  1 

generates  the  delay  time  t  ,  .   This  output  synchronizes  the 
generation  of  the  next  pulse  or  oscillation  at  time  T1  which 
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reach  their  receivers  after  time  intervals  t'  ,  and  t'  „ 

wl        w2 

seconds.   Therefore,  the  pulse  spacing  is  t    and  t    and 

the  total  time  from  the  emitting  of  the  first  pulse,  T,  to 

the  receiving  of  the  second  pulse,  T1.  by  R  .  is  t  .  +  t  , 

w2      wl     Dl 

+  t1  ,. 

w2 

In  terms  of  the  speed  of  sound  in  still  air,  C,  and 
the  separation,  % ,  of  the  transmitter  and  receiver, 


21 
wl     Dl      w2    C      Dl 


(15) 


As  previously  discussed,  the  temperature  measurement  is 

sampled  at  half  the  rate  of  the  wind  measurement  (220  times/ 

sec  as  compared  with  440  times/sec).   Therefore,  the  time, 

t  =  t  ,  +  t'  o  is  related  to  temperature.   Since  Atn  changes 
wl      w2  c  8 

very  slowly  with  temperature,  a  time  t  -  generated  by 

another  time  delay  circuit  is  subtracted  from  t  +  t  1 .   The 

time  difference,  AtQ  =  t  +  fc  _  -  t^-,  is  then  supplied  to  a 

0  Dl     D  Z 

temperature  analog  card  which  changes  it  into  a  Atq  (the 

mean  temperature  found  on  the  8  dial  of  the  temperature 

unit.   The  dial  is  operated  manually).   The  time  difference 

is  also  changed  into  At'  ,  the  fluctuating  component,  which 

is  changed  into  a  voltage  that  appears  on  the  0  y    meter 

(Figure  14) . 

The  time  delay  t  _  is  adjustable  (longer  at  lower 

temperatures)  since  it  is  preferable  to  keep  Atfi  small.   If 

At0  is  made  too  small,  however,  a  sluggish  limiting  of  the 
U 

temperature  signal  will  occur  —  the  effect  showing  as  a 
rounding  to  flattening  of  the  high  temperature  peaks. 
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Figure  13.   Timing  Sequence  Schematic.  Kaijo  Denki  (1969) 
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Figure  14.   Temperature  and  Time  Relationship 
Kaijo  Denki  (1969) 
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III.   DATA  ANALYSIS 

A.   GENERAL 

Turbulent  processes  appear  as  erratic  changes  in  tem- 
perature and  in  direction  and  speed  of  the  wind.   A  know- 
ledge of  the  nature  of  turbulent  processes  in  the  atmosphere 
(fluctuations  of  such  parameters  as  wind  velocity,  tempera- 
ture and  moisture  for  example)  can  lead  to  a  better  under- 
standing of  heat,  energy  and  momentum  transfers  which  occur 
over  the  ocean. 

The  mechanics  of  turbulence  are  so  complex  that  it  is 
almost  traditional  now  to  use  statistical  methods  through 
which  an  understanding  of  the  mean  characteristics  can  be 
obtained  if  not  in  detail.   One  statistical  method  which 
has  been  used  with  some  considerable  success  is  spectral 
analysis.   From  spectral  analysis  the  distribution  of  energy 
of  a  fluctuating,  turbulent  process  (e.g.  wind  velocity, 
temperature)  can  be  plotted  as  a  function  of  frequency. 

Analog  measurements  were  recorded  of  wind  velocity  and 
temperature  fluctuations  from  the  sonic  anemomet er/ thermo- 
meter, hot  wire  anemometer  and  a  one  mil  thermocouple  then 
digitized  and  then  spectral  analyzed  by  a  Fast  Fourier 
Transform  method.   A  general  discussion  of  procedures  used 
in  the  analog-to-digital  conversion  and  spectral  analysis 
will  be  given  in  the  next  sections.   More  detailed  discus- 
sion of  procedures  and  theory  appear  in  McKendrick  (19  71) 
Jones  (1971)  and  Wilson  et  al  (1969). 
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B.   KOLMOGOROV'S  HYPOTHESIS 

1 .   General 

An  important  consideration  in  turbulence  theory  is 

whether  the  turbulence  is  anisotropic  or  isotropic.   Aniso- 

tropy  exists  for  large  scale  eddies  and  isotropy  is  thought 

to  exist  for  small  scale  eddies.   Isotropy  was  described  by 

Lumley  and  Panofsky  (1964)  as  follows: 

"An  isotropic  field  is  one  in  which  there  is  no 
preferred  direction.   Thus  all  average  functions  de- 
scribing the  statistics  of  such  a  field  must  remain 
unchanged  regardless  of  a  rotation  (or  reflection)  of 
coordinate  system.   For  example,  the  spectrum  of  a 
function,  (J)  (k)  ,  homogeneous  in  three  spatial  directions, 
which  is  also  isotropic,  cannot  be  a  function  of  the 
separate  components  of  the  wave  number  vector,  for  a 
rotation  to  a  new  coordinate  system  would  change  those. 
It  can  be  a  function  only  of  the  length  of  the  vector, 
for  that  is  the  only  quantity  characterizing  'k'  which 
does  remain  unchanged  under  a  rotation." 

Kolmogorov  postulated  that  at  high  wave  numbers 
turbulence  should  be  isotropic  even  though  it  is  generated 
from  and  embedded  in  anisotropic  eddies  of  low  wave  number-. 
If  the  Reynolds  number  (ratio  of  inertial  forces  to  viscou 
forces  in  the  flow)  is  high  enough,  turbulence  will  adjust 
until  a  statistically  steady  state,  independent  of  initial 
conditions,  is  attained.   Turbulence  is  adjusted  through 
inertial  transfer  (forces  due  to  fluctuating  velocity  gradi- 
ents) and  viscous  dissipation.   This  steady  state  consists 
of  energy  transfer  by  inertial  forces  into  the  high  wave 
number  range  that  is  equal  to  energy  lost  through  viscous 
dissipation.   This  predicted  behavior  of  velocity  spectra 
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was  utilized  in  interpretation  and  comparison  of  the  sonic 
anemomeme t er ,  hot  wire  anemometer  and  one  mil  thermocouple 
data.  Prediction  by  Kolmogorov's  hypotheses  are  discussed 
in  more  detail  in  the  following  paragraphs. 

2 .   Kolmogorov's  First  Hypothesis:  Equilibrium  Range 

At  high  Reynolds  numbers  "small  scale"  eddies  are 
isotropic  and  in  universal  statistical  equilibrium.   Their 
properties  depend  on  dissipation  (e)  and  on  kinematic 
viscosity  (v)  and  not  directly  on  time.   Motion  associated 
with  the  equilibrium  range  is  uniquely  determined  by  £ 
and  V  . 

Kolmogorov  concluded  that  the  energy  spectrum 
[E(k)],  in  units  of  energy  per  unit  wave  number,  is 

,   5,1/4     *  ^kc; 
(ev  )  S 


where 


f  (-r — )  =  undesignated  function 


S  . 


■  <7->lM 


~    characteristic  wave  number 


3 .   Kolomogor ov ' s  Second  Hypothesis:  Inertial  Subrange 
If  the  Reynolds  number  is  large  enough,  the  energy 
containing  eddy  wave  number  region  and  the  "dissipation" 
wave  number  region  should  be  well  separated  in  wave  number 
space.   There  should  exist  a  range  of  wave  numbers  (eddy 
sizes)  on  the  low  "k"  side  of  the  dissipation  range  (but  at 
higher  wave  number  than  energy  containing  wave  numbers)  in 
which  the  universal  statistical  steady  state  holds 


3  6 


["x~r    E(k)  =  0]  but  where  there  is  negligible  dissipation  of 

o  t 

energy,  [Figure  15  from  Lumley  and  Panofsky  (1964)].  This 
range  of  wave  numbers  where  the  dominant  feature  is  energy 
flux  by  inertial  forces  is  called  the  inertial  subrange. 

Since  viscosity  is  not  important  in  this  wave  number 
range,  statistical  properties  are  determined  only  by  dissi- 
pation rate  (e)  and  wave  number  (k) .   Kolmogorov  stated  that 
energy  spectrum  should  have  the  following  functional  rela- 
tion 

E(k)  -  C  e2/3  k"3/3 

£  =  dissipation  of  turbulent 
kinetic  energy 

k  =  wave  number 

C  =  constant  of  proportionality 
Since  frequency  is  directly  proportional  to  wave 
number,  a  log-log  plot  (base  10)  of  energy  spectra  versus 
frequency  would  yield  a  curve  with  a  -5/3  slope  if  the  spec- 
tral estimates  are  within  the  frequency  band  of  the  inertial 
sub  range . 

C.   ANALOG  TO  DIGITAL  PROCEDURES 
1 .   General 

Analog  data  recorded  in  the  field  were  digitized 
through  use  of  a  Hybrid  computer  system.    The  Hybrid  system 
consists  of  an  analog  computer,  C0MC0R  Ci-5000,  electron- 
ically connected  to  a  digital  computer,  XDS-9300.   A  program 


This  system  is  located  in  Spanagel  Hall,  Naval  Postgraduate 
School,  Monterey,  California. 
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Figure  15.   Schematic  Drawing  of  Spectral  Energy  Transfer 
in  Free  Convection.   Lumley  and  Panof sky  (1964) 
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exists  for  digitizing,  using  the  XDS-9300  and  the  Ci-5000 
analog  computer.   An  analog  patchboard  and  a  logic  board, 
which  are  plugged  into  the  Ci-5000  were  also  available. 
Digitization  was  performed  in  the  octal  number  base  onto 
seven  track  magnetic  tapes. 

2  .   Considerations  Prior  to  Digitization 

Considerations  taken  into  account  prior  to  digiti- 
zation were  noise,  D.C.  level,  analog  signal  strength  and 
total  analog  data  record  length,  T  (corresponding  to  time). 
All  of  these  were  functions  of  collection  technique. 

Noise  that  contaminates  the  desired  signal  is 
usually  of  high  frequency  and  normally  can  be  low  pass  fil- 
tered during  the  data  recording  procedures.   This  is  the 
preferred  method.   However,  filtering  can  be  done  in  the 
Hybrid  computer  facility  just  prior  to  digitizing.   In  this 
study,  the  analog  data  were  filtered  at  the  Hybrid  computer 
Data  from  the  four  sonic  anemomet er/ thermometer  channels 
were  filtered  at  31.25  Hz  and  the  hot  wire  anemometer  data 
and  one  mil  thermocouple  data  were  filtered  at  100.0  Hz. 
(Filter  frequencies  should  equal  one-half  of  the  digital 
sampling  rate .  ) 

The  D.C.  level  was  removed  from  all  data  during 
digitization.   If  it  were  necessary  that  the  D.C.  level  be 
removed  from  the  data  during  data  recording,  a  band  pass 
filter  with  a  low  frequency  cut  off  of  less  than  one  hertz 
could  have  been  used.   This  would  preserve  some  of  the  low 
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frequencies  while  eliminating  the  high  frequency  noise  and 
the  D.C.  level. 

All  analog  signals  were  amplified  in  order  to  take 
full  advantage  of  the  range  of  the  Hybrid  computer  system. 
The  Hybrid  system  is  able  to  digitize  signals  of  +  100  volts. 
Channels  A  and  B  of  the  sonic  anemometer/ thermometer  and 
the  hot  wire  anemometer  were  amplified  by  50;  Channels  W 
and  T  of  the  sonic  anemometer/ thermometer  and  the  thermo- 
couple were  amplified  by  100.   Amplification  factors  were 
taken  into  account,  in  this  study,  in  the  Fast  Fourier  Trans- 
form program. 

Consideration  had  to  be  given  to  obtaining  an  opti- 
mum analog  data  record  length,  T,  during  the  field  experi- 
ments.  In  optimizing  this  total  record  length  for  the 
frequencies  of  interest,  the  relationship  between  T  and  the 
shape  of  the  convolved  wave  forms  (Fourier  Transform)  must 
be  considered.   Jones  (1971)  described  fully  the  effects  of 
record  lengths  on  the  shape  of  the  convolved  wave  forms. 
In  general  there  is  some  optimum  record  length  where  the 
difference  between  the  convolved  spectrum  and  true  spectrum 
is  very  small.   Although  Jones  (1971)  did  not  find  an  opti- 
mum record  length,  he  found  that  increasing  a  record  from 
140  seconds  to  360  seconds  had  no  significant  effect  on 
spectral  shape  between  10  Hz  and  1000  Hz.   Data  collected  in 
this  study  had  a  record  length  of  30  minutes. 

Low  frequency  resolution  is  also  tied  directly  to 
the  length  of  the  analog  record.   According  to  Rayleigh's 


40 


Criterion,  it  is  hard  to  distinguish  low  frequencies  (f)  in 
the  interval  0  <  f  <  1/2T  (where  T  =  analog  record  length 
in  seconds).   When  dealing  with  a  finite  record  the  original 
signal  is  truncated.   A  long  record  narrows  the  frequency 
resolution.   Record  lengths  of  30  minutes  jave  a  low  fre- 
quency resolution  of  approximately  0.0003  Hz  although  sta- 
tistical reliability  would  not  be  expected  until  0.003  Hz. 
Another  factor  which  had  to  be  considered  prior  to 
digitizing  was  that  of  aliasing.   Aliasing  is  dependent  on 
the  folding  (Nyquist)  frequency  which  is  1/2  the  sampling 
frequency.   In  order  to  prevent  the"folding  back"  of  energy 
from  higher  frequencies  into  the  lower  frequencies  of 
interest,  a  sampling  rate  greater  than  twice  the  highest 
frequency  desired  must  be  chosen.   The  folded  back  energy 
will  be  centered  around  the  Nyquist  frequency.   Sonic  anemo- 
meter/thermometer data  from  Experiment  I  was  run  at  two 
different  sampling  rates.  Sampling  rates  of  390.5  Hz  and 
62.5  Hz  gave  aliasing  frequencies  of  195.25  Hz  and  31.25  Hr 
respectively.   The  hot  wire  anemometer  data  and  the  thermo 
couple  data  from  Experiment  II  were  sampled  at  200.0  Hz  which 
gave  an  aliasing  frequency  of  100.0  Hz.   Experiment  II  sonic 
anemometer/ thermometer  data  were  sampled  at  62.5  Hz  only. 

D.   CONVERSION  OF  SEVEN  TRACK  TO  NINE  TRACK 

Digitization  of  analog  data  on  the  above  system  produced 

a  seven  track  tape  in  the  octal  number  c?se.  To  utilize 

this  data  in  the  IBM-360/67  compucer  system,  digitized  data 
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had  to  be  converted  to  a  nine  track  tape  in  the  hexadecimal 
number  base.   This  conversion  process  was  done  using  an 
available  program,  CONVERT.   Complete  listings  and  Job  Con- 
trol Language  (JCL)  card  setups  for  the  CONVERT  program  can 
be  found  in  theses  by  McKendrick  (1971)  and  Jones  (1971), 
and  in  a  report  by  Raney  (1971). 

E.   FOURIER  COEFFICIENT  COMPUTATIONS 

Fourier  coefficients  were  obtained  by  using  the  FTOR 
program,  one  of  two  UBC  time  series  programs  currently  on 
file  in  the  IBM-360/67  subroutine  library  at  the  Naval  Post- 
graduate School. 

Transformation  of  digitized  signal  data  from  the  time 
domain  into  the  frequency  domain  was  done  using  a  finite 
form  of  the  Fourier  Transform.   Data  were  read  off  the  con- 
verted nine  track  tape,  a  record  at  a  time.   The  Fourier 
coefficients  were  computed  and  were  then  stored  on  a  sepa- 
rate tape  . 

The  longer  the  record,  the  closer  together  in  frequency 
will  be  the  spectral  estimates  because  the  value  of  the 
frequency  is  given  by  (n/T)  .   That  is, the  discrete  frequency 
bandwidth  is  inversely  proportional  to  the  record  length. 

Due  to  the  prohibitive  amount  of  computer  time  used  to 
compute  coefficients,  the  direct  finite  Fourier  Transform 
of  a  series  of  data  points  was  not  used.  The  Fast  Fourier 
Transform  of  Cooley  and  Tukey  (1965)  was  used  and  it  is  a 
variation  of  the  direct  method.   The  Fast  Fourier  Transform 
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is  incorporated  into  the  FTOR  program  resulting  in  consider- 
ably less  computer  time  for  Fourier  coefficient  computation. 
This  method  yields  an  improvement  in  excess  of  five  orders 
of  magnitude  in  computational  time  (McKendrick  1971)  over 

the  direct  method. 

13 

The  FTOR  program  transforms  a  maximum  of  2   (8192)  data 

points.   This  maximum  is  set  by  dimension  limitations  in  the 
program.   Turbulence  data  normally  includes  many  more  data 
points  than  this  so  the  data  points  are  read  in  as  records 

consisting  of  integral  powers  of  two  (up  to  a  maximum  of 

13 
2   )  and  these  records  are  transformed  consecutively. 

F.   UBCSCOR  PROGRAM 

Once  the  Fourier  coefficients  have  been  computed,  the 

second  UBC  program,  SCOR,  converts  them  to  power  spectral 

densit  ies . 

Spectral  values  E.,(n)  are  computed  from  the  complex 

3  k 

Fourier  coefficients  by  forming  the  product  of  the  complex 

coefficient  [A  (n)  =  a . (n)  -  i  b.(n)]  with  a  conjugate 

* 

[A   (n)  =  a,  (n)  +  i  b,  (n) ]  where  j  and  k  could  be  the  same 

or  different . 

Spectral  values  were  obtained  by  the  following  computa- 
t  ions 


Spectrum  density:   j  =  k, 


E. . (n)  = 
J  J 


A  (n) »A  (n)        .      ?   t 

J  9A— J =  (a    +  b  Z)  y 

2£n  n      n    2 


where  An  =  — 
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Cross  spectrum :  j  ^  k , 

A.(n) -A*(n) 
E.,  (n)  =  -L 


2An 

=  ■£-[  (a  c  +b  d  )-i(b  c  -a  d  )  ] 
2    nn   nn      nn   n  n 

The  real  and  imaginary  parts  of  the  cross  spectrum  are 

the  cospectrum  and  quadrature  spectrum  respectively  and 

were  denoted  as;     Cospectrum : 

E..  (n)  =  (a  c   +  b  d  )  \ 
jk         nn     n  n   2 

Quadrature  spectrum: 


E.,  (n)  =  (b  c   -  a  d  )  -i- 
j  k         n  n     n  n   2 


The  spectra  and  cospectra  are  computed  for  each  fre- 
quency in  successive  records  and  then  these  records  are 
averaged  consecutively.   The  program  also  divides  the  spec- 
tral densities  into  32  bandwidths  and  averages  them  across 
the  bandwidths  as  well  as  along  the  records  to  obtain  a 
mean  and  a  variance  of  the  quantities. 

In  the  SCOR  program  one  can  specify  the  number  of 
channels  to  be  analyzed,  the  number  of  records  in  a  file, 
the  start  point  for  analysis,  selection  of  spectra  and/or 
cospectra,  bandwidth  desired  (logarithmic  or  linear),  the 
type  of  graphic  plot  desired  and  for  this  study,  machine-set 
coordinates  for  the  graphic  output.   Pro  gramme r-se t -co- 
ordinates also  can  be  used.   The  program  computes,  the 
spectra  and/or  cospectra  for  each  record  in  sequence  until 
all  records  requested  have  been  analyzed  or  the  "end-of- 
file"  mark  is  reached  on  the  tape. 
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IV  .        FIELD    EXPERIMENTS 

A.  PRELIMINARIES 

The  sonic  anemometer/thermometer  used  in  the  field  ex- 
periments had  been  in  storage  for  a  period  of  two  years. 
The  first  consideration  was  to  thoroughly  check  the  sonic 
anemometer/ thermometer  electronically.   It  was  checked  out 
in  accordance  with  manufacturer's  instructions  [Kaijo  Denki 
(196 9 ), Append ix  D  and  E  which  included  voltage  and  wave  form 
checks  at  certain  critical  points].   A  large  number  of  tran- 
sistors in  both  the  power  supply  and  junction  box  had  to 
be  replaced . 

Once  the  sonic  anemometer / thermomet er  was  checked  elec- 
tronically, wind  tunnel  tests  were  performed.   Simple, 
operational  tests  in  a  controlled  environment  were  made  to 
insure  that  the  sonic  anemomet er /thermomet er  was  measuring 
wind  speed  and  temperature  correctly.   Once  satisfactory  wind 
tunnel  measurements  were  made,  field  experiments  were 
conducted . 

B.  OBJECTIVE  OF  EXPERIMENTS 

The  objective  of  the  field  experiments  was  to  gather  com- 
parative turbulence  data  from  the  sonic  anemomet er/ thermo- 
meter, the  hot  wire  anemometer  and  a  one  mil  thermocouple. 
All  three  sets  of  data  were  then  similarly  spectrally 
analyz  ed . 
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C.   EXPERIMENTS 

1 .  Installation  of  Instruments 

The  procedures  described  in  this  section  applies  to 
both  experiments.   All  equipment  was  installed  at  the  Naval 
Postgraduate  School  beach  lab , indicated  in  Figure  16  and 
Figure  17.   The  sonic  probe,  hot  wire  probe  and  one  mil  ther- 
mocouple were  mounted  on  a  pole  at  a  height  of  three  meters 
(Figure  18)  and  positioned  so  that  all  sensors  were  facing 
directly  into  the  wind  flow  which  was  blowing  off  the  water. 
During  the  experiments,  all  probes  were  kept  into  the  wind 
such  that  the  velocity  measured  on  channels  A  and  B  were 
equal.   A  hand  held  wind  vane  was  used  to  check  on  the  wind 
direction  in  order  to  maintain  proper  probe  orientation. 

All  recording  dev ices , power  supplies  and  amplifiers 
were  installed  inside  the  beach  lab  (Figure  19) .   Analog 
data  were  recorded  on  a  SANGAMO  FM  tape  recorder  model  3562 
at  a  speed  of  1-7/8  inches  per  second.   This  recorder  uses 
one  inch  magnetic  tape  and  has  14  tracks  available  for  data 
and  two  edge  tracks  available  for  voice  recording. 

The  hot  wire  anemometer  used  for  velocity  compari- 
sons was  a  TSI  model  1050.   A  copper  constantan,  unsheathed, 
fine  gauge  Omego  Engineering  Company  ( SSTSSP-02 OE-6 )  thermo- 
couple was  used  for  temperature  measurements. 

2 .  Experiment  I  -  14  December  1972 

Wind  speed  and  temperature  measurements  were  made 
during  the  late  afternoon  on  14  December  1972.   Visibility 
was  unlimited.   Winds  were  from  the  west  to  northwest,  off 
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Figure  16.   Arrangement  of  Probes  for  Experiments  I  and  II 
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Figure  17.   Closeup  of  Probe  Arrangement 
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Figure  18.   Probe  Stand  Orientation 
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Figure  19.   Inside  Installation  of  Equipment 
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the  ocean,  at  an  average  speed  of  2-4  m/sec.   Mean  air  tem- 
perature for  the  period  (1430-1700  PST)  was  52°F. 

Three  separate  30  minute  analog  data  records  were 
recorded  (1500  to  1530,  1546  to  1616  and  1620  to  1650  PST). 
Data  on  all  channels  were  recorded  without  amplification. 
During  the  measurements,  spikes  were  observed  on  all  channels 
The  spikes  were  found,  after  this  experiment,  to  have  been 
caused  by  faulty  transistors  in  the  junction  box.   Examples 
of  these  spikes  appear  in  the  strip  chart  recording  in 
Figure  20. 

Hot  wire  anemometer  and  thermocouple  data  were  not 
recorded  because  the  record  amplifier  on  the  hot  wire  channel 
of  the  tape  recorder  failed  and  the  thermocouple  malfunc- 
tioned . 

The  three  separate  analog  data  records  from  the 
sonic  anemomete r / thermome ter  data  were  digitized  and  spectral 
analyzed . 

3.   Experiment  II  -  27  February  1973 

This  experiment  was  similar  to  experiment  I  except 
that  hot  wire  anemometer  data  and  thermocouple  data  were 
successfully  recorded.   Channel  B  of  the  sonic  anemometer/ 
thermometer  was  not  operating  because  of  a  short  in  the  con- 
nector cable  on  that  channel.   This  malfunction  did  not 
affect  the  other  channels. 

Again,  wind  speed  and  temperature  data  were  obtained 
during  three  separate  30  minute  periods  (1030  to  1100,  1345 
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Figure  20.   Strip  Chart  Recording  of  the  Sonic  A  and  W 
Channels  Indicating  Spikes. 
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to  1415  and  1530  to  1600  PST).   It  was  an  overcast  sky. 
Visibility  was  approximately  six  miles.   Winds  were  from  the 
west,  off  the  ocean,  with  an  average  speed  of  4-8  m/sec. 
Mean  air  temperature  ranged  from  58°F  in  the  late  morning 
period  to  60°F  in  the  early  afternoon  period. 

During  the  recording  periods,  all  but  one  of  the 
data  channels  were  amplified  using  Preston  8300-XWB  ampli- 
fiers.  The  amplification  factors  were  two  for  the  sonic 
W  and  T  channels,  one  for  the  sonic  A  channel,  twenty  for 
the  hot  wire  channel  and  one  thousand  for  the  thermocouple- 
channel.   The  hot  wire  data  were  filtered  before  recording 
using  a  Krohn-Hite  model  3750  bandpass  filter  low  set  at 
0.2  Hz  and  high  set  at  200,0  Hz.   The  three  separate  30 
minute  analog  data  records  were  digitized  and  spectrally 
analyzed . 
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V.   RESULTS  AND  CONCLUSIONS 

A.  GENERAL  PERFORMANCE  OF  SONIC  ANEMOMETER/THERMOMETER 
In  general  the  sonic  anemometer /thermometer  satisfac- 
torily measured  the  wind  and  temperature  fluctuations  during 
the  field  experiments,  but  it  did  not  operate  at  its  peak 
potential  because  of  various  electronic  problems.   For 
example,  spiking  occurred  on  all  channels  during  both  ex- 
periments I  and  II  and  a  short  occurred  in  channel  B  during 
experiment  II.   These  could  have  been  a  result  of  the  instru- 
ment being  idle  for  two  years. 

Spiking  proved  to  be  a  continuous  problem  during  the 
field  experiments  even  though  numerous  transistors  in  the 
system  had  been  replaced. 

B.  SPECTRA  PLOTS 
1 .   General 

Sonic  anemometer/ thermometer  data  recorded  on  14 
December  1972  (experiment  I)  were  digitized  at  a  sampling 
rate  of  390.5  Hz  (aliasing  frequency  192.25  Hz)  and  also 
at  62.5  Hz  (aliasing  frequency  31.25  Hz. 

No  editing  or  filtering  was  performed  on  any  of  the 
channels  during  digitization  at  the  higher  sampling  rate. 
The  higher  sampling  rate  (390.5  Hz)  enabled  the  examination 
of  a  frequency  band  from  0  Hz  to  approximately  100  Hz. 

The  analog  data  tape  was  edited  prior  to  digitizing 
at  the  lower  sampling  rate.   This  led  to  the  elimination  of 
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the  first  time  period  which  had  the  heaviest  spiking.   The 
remaining  two  time  periods  were  digitized  at  62.5  Hz  and 
filtered  at  31.25  Hz.  The  lower  sampling  rate  shifted  the 
frequency  band  examined  (approximately  0.5  Hz  to  40  Hz) 
towards  the  lower  end  of  the  frequency  spectrum.  The  spec- 
tral results  indicate  both  the  energy  containing  frequency 
band  as  well  as  the  lower  frequency  er.d  of  the  inertial  sub- 
range . 

Hot  wire  data  were  not  taken  during  experiment  I 
because  of  a  recording  amplifier  failure.   Two  samples  of 
hot  wire  data  taken  on  a  day  that  had  similar  synoptic 
characteristics  (air  temperature  and  wind)  as  the  day  of 
experiment  I  are  included  for  comparison  purposes.   They 
are  considered  representative  samples  of  data  which  would 
have  been  taken  by  a  hot  wire  anemometer  system  on  that  day. 
These  results  represent  a  frequency  band  of  from  0.01  Hz 
to  approximately  100  Hz.   The  data  for  both  of  these  results 
were  digitized  at  a  rate  of  200.0  Hz  (aliasing  frequency 
100.0  Hz)  and  filtered  at  100.0  Hz. 

On  28  February  1973  (experiment  II),  sonic  anemo- 
meter/thermometer, hot  wire  anemometer  and  one  mil  thermo- 
couple data  were  collected  during  the  three  30  minute  time 
periods.   The  sonic  data  were  digitized  at  a  sampling  rate 
of  62.5  Hz  and  filtered  at  32.25  Hz.   The  hot  wire  and 
thermocouple  data  were  digitized  at  a  sampling  rate  of 
200.0  Hz  and  filtered  at  100.0  Hz.   Sonic  channel  B  was 
inoperat  ive . 
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2  .   Sonic  Anemometer/Thermometer  Data  Sampled  at 
390.5  Hz  -  14  December  1972  (Experiment  I) 

Figure  21  to  Figure  32  are  for  the  three  30  minute 
periods  from  1500  to  1700.   A  -5/3  slope  exists  on  the  log- 
log  plot  in  the  frequency  band  from  1.0  Hz  to  10.0  Hz  for 
all  channels.   Spectra  from  channels  A,  B,  W  have  -5/3 
slopes  which  continue  beyond  10.0  Hz  but  the  T  spectra 
slopes  drop  off  more  rapidly  than  -5/3  before  10.0  Hz. 

3 .  Sonic  Anemometer/Thermometer  Data  Sampled  at 
62.5  Hz  -  14  December  1972  (Experiment  I) 

Figure  33  to  Figure  40  are  from  the  two  edited 

periods  of  data  which  were  filtered  and  then  sampled  at  the 

lower  rate.    Spectra  from  channels  A  and  B  in  both  periods 

have  definite  -5/3  slopes  within  the  1.0  Hz  to  10.0  Hz  band 

with  some  departure  around  8.0  Hz.   For  both  periods,  W 

channel  spectra  have  -5/3  slopes  up  to  6.0  Hz  with  rapid 

deterioration  beyond.   For  the  second  period,  the  T  spectrum 

has  a  drop  off  from  a  -5/3  slope  at  approximately  4.0  Hz  but 

for  the  third  period,  the  T  spectra  has  a  definite  -5/3 

slope  throughout  the  1.0  Hz  to  10.0  Hz  frequency  band. 

4 .  Hot  Wire  Anemometer  Data  Sampled  at  200.0  Hz 
-  10  January  1973 

Figure  41  and  Figure  42  are  from  hot  wire  anemometer 

data  taken  on  10  January  1973.   Two  periods  are  discussed 

and  each  period  was  approximately  50  minutes  in  length.  A 

distinct  -5/3  slope  is  found  in  the  frequency  band  from 

1.0  Hz  to  approximately  40.0  Hz  for  both  periods  of  data. 
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Figure  21.   Sonic  "A"  Channel  Data  Observed  14  December 
1972  1st  Period  (Experiment  I)  High  Sample  Race. 
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Figure  22.   Sonic  "B"  Channel  Data  Observed  14  December 
1972  1st  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  23.   Sonic  "W"  Channel  Data  Observed  14  December 
1972  1st  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  24.   Sonic  "T"  Channel  Data  Observed  14  December 
1972  1st  Period  (Experiment  I)  High  Sample  Rate. 


60 


(VJ 


q: 

LlJ 

LEI 
S   o 


10HZ 


ID 
CE 
h- 
U 
LU 

o 


100HZ 


-0.  <1 


LOGIC   FREQUENCY 


i.  6 


2.  4 


Figure  25.   Sonic  "A"  Channel  Data  Observed  14  December 
1972   2nd  Period  (Experiment  I)  High  Sample  Rate. 


61 


r\i 


en 

LlJ 

x 
s 


(\l 


o 

LjJ 

O 
•H 
Li 

_l 

a 

19 


-3- =3 


10HZ 


100HZ 


L.aC.10   FREQUENT 


Figure    26.       Sonic    "B"    Channel    Data    Observed    14    December 
1972    2nd    Period     (Experiment    I)    High    Sample    Rate. 
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Figure  27.   Sonic  "W"  Channel  Data  Observed  14  December 
1972  2nd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  28.   Sonic  "T"  Channel  Data  Observed  14  December 
1972  2nd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  29.   Sonic  "A"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  30.   Sonic  "B"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  31.   Sonic  "W"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  32.   Sonic  "T"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  High  Sample  Rate. 
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Figure  33.   Sonic  "A"  Channel  Data  Observed  14  December 
1972  2nd  Period  (Experiment  I)  Low  Sample  Rate. 
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"B"  Channel  Data  Observed  14  December 


Figure  34 .   Sonic 

1972  2nd  Period  (Experiment  I)  Low  Sample  Rate. 
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Figure    35.       Sonic    "W"    Channel    Data    Observed    14    December 
1972    2nd    Period     (Experiment    I)    Low    Sample    Rate. 
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Figure  36.   Sonic  "T"  Channel  Data  Observed  14  December 
1972  2nd  Period  (Experiment  I)  Low  Sample  Rate. 
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Figure  37.   Sonic  "A"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  Low  Sample  Rate. 
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Figure  38.   Sonic  "B"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  Low  Sample  Rate. 


74 


h- 

cr 

UJ 

x 
\ 

o 


cc 
\— 

LU 

D_  o 

t 

D 
—I 
LD 


1.0  HZ 


100  HZ 


-j.  5  -a.  9  0.  a  :.  3 

LGG10   FREQUENCE 


1.6 


Figure  39.   Sonic  "W"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  Low  Sample  Rate. 
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Figure  40.   Sonic  "T"  Channel  Data  Observed  14  December 
1972  3rd  Period  (Experiment  I)  Low  Sample  Rate. 
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Figure  41.   Hot  Wire  Data  Observed  10  January  1973 
1st  Period . 
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Figure  42.   Hot  Wire  Data  Observed  10  January  1973 
2nd  Period  . 
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5 .  Sonic  Anemome t er /Thermome t er  Data  Sampled  at 
62.5  Hz  -  28  February  1973  (Experiment  II) 

Figure  43  to  Figure  51  are  spectra  from  the  sonic 

anemometer/ thermometer  data  recorded  on  28  February  1973. 

Channel  B  was  inoperative  during  the  period.   Spectra  from 

channels  A  and  W  show  a  definite  -5/3  slope  in  the  1.0  Hz 

to  10.0  Hz  frequency  band.   A  -5/3  slope  in  a  very  narrow 

band  of  frequencies  of  approximately  1.0  Hz  to  4.0  Hz 

appears  in  the  T  spectra. 

6 .  Hot  Wire  Anemometer  Data  Sampled  at  200.0  Hz 
-  28  February  1973  (Experiment  II) 

Figure  52  and  Figure  53  are  hot  wire  spectra  from 

the  second  and  third  periods  recorded  on  28  February  1973. 

The  first  period  of  data  was  lost  due  to  a  faulty  magnetic 

tape.   Spectra  from  both  periods  have  definite  -5/3  slopes 

in  the  region  1.0  Hz  to  approximately  91.0  Hz. 

7 .  Thermocouple  Data  Sampled  at  200.0  Hz  - 
28  February  1973  (Experiment  II) 

Figure  54  and  Figure  55  are  thermocouple  spectra 

from  the  second  and  third  periods  recorded  on  28  February 

1973.  The  first  period  of  data  in  this  case  was  also  lost 

due  to  the  faulty  magnetic  tape.   From  the  second  period 

plot  the  -5/3  slope  starts  at  approximately  4.4  Hz  and  goes 

just  beyond  10.0  Hz.   From  the  third  period  plot,  a 

definite  -5/3  slope  exists  in  the  range  1.0  Hz  to  10.0  Hz. 

C.   CONCLUSIONS 

The  sonic  anemome te r/ thermome ter  system  proved  to  be  a 
reliable,  rugged  and  easily  installed  instrument  for  the 
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Figure  43.   Sonic  "A"  Channel  Data  Observed  28  February 
1973  1st  Period  (Experiment  II). 
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Figure  44.   Sonic  "W"  Channel  Data  Observed  28  February 
1973  1st  Period  (Experiment  II). 
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Figure  45.   Sonic  "T"  Channel  Data  Observed  28  February 
1973  1st  Period  (Experiment  II). 
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Figure  46.   Sonic  "A"  Channel  Data  Observed  28  February 
1973  2nd  Period  (Experiment  II)  . 
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Figure  47.   Sonic  "W"  Channel  Data  Observed  28  February 
1973  2nd  Period  (Experiment  II)  . 
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Figure  48.   Sonic  "T"  Channel  Data  Observed  28  February 
1973  2nd  Period  (Experiment  II). 
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Figure  49.   Sonic  "A"  Channel  Data  Observed  28  February 
1973  3rd  Period  (Experiment  II) . 
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Figure    50.        Sonic    "W"    Channel    Data    Observed    28    February 
1973    3rd    Period     (Experiment    II) . 
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Figure  51.   Sonic  "I"  Channel  Data  Observed  28  February 
1973  3rd  Period  (Experiment  II). 
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Figure  52.   Hot  Wire  Data  Observed  28  February  1973 
2nd  Period  (Experiment  II). 
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Figure  53.   Hot  Wire  Data  Observed  28  February  1973 
3rd  Period  (Experiment  II). 
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Figure  54.   Thermocouple  Data  Observed  28  February  1973 
2nd  Period  (Experiment  II). 
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Figure  55.   Thermocouple  Data  Observed  28  February  1973 
3rd  Period  (Experiment  II). 
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measurement  of  the  three  components  of  the  wind  and  also  air 
temperature.   Its  frequency  response  was  found  to  be  from 
0.0  Hz  to  approximately  12.0  Hz. 

The  sonic  anemometer/ thermometer  measurements  did 
not  appear  to  be  affected  by  weather  elements  and  because  it 
uses  the  speed  of  sound  in  air  it  has  inherent  calibration. 
The  response  of  the  sonic  anemome t er / the rmome t er  is  linear 
and  so  is  its  sensitivity.   It  is  thus  able  to  be  used  to 
describe  vertical  and  horizontal  transfers  of  energy,  heat 
and  momentum. 

The  -5/3  slope  of  Kolmogorov's  inertial  subrange 
was  easily  recognizable  on  spectral  plots  of  turbulence  data 
measured  by  the  sonic  anemometer/ thermome ter .   Thus  the  com- 
putation of  the  rate  of  dissipation  of  kinetic  energy  con- 
stant can  be  made  which  in  turn  means  that  heat  and  momentum 
transfer  computations  can  be  made. 
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APPENDIX  A 
Definition  of  Terms  Used  in  Analog-Digital  Data  Processing 

The  terms  used  in  the  analog- to-digital  process  are  de- 
fined here  based  on  current  usage  in  the  Hybrid  computer 
center  and  based  on  that  utilized  by  the  three  time  series 
programs  for  the  spectral  analysis  of  digitized  data.   There 
are  some  contradictions  of  terms  between  McKendrick  (1971) 
and  Jones  (1971)  and  between  what  is  currently  used. 

The  word  "record"  is  often  used  in  both  analog  and  digi- 
tal discussions  as  meaning  a  finite  length  of  data  in 
seconds.   In  the  analog  sense,  "record"  means  the  entire 
length  of  data  which  is  to  be  digitized.   In  the  digital 
sense,  "record"  means  a  short  length  of  data  based  on  the 
sampling  rate  and  the  number  of  samples  per  record  selected 
by  the  programmer.   In  McKendrick  (1971)  and  Jones  (1971), 
digital  "record"  is  referred  to  as  "block".   The  word 
"block"  is  not  used  in  this  thesis  as  "record"  is  what  is 
utilized  in  all  three  programs  (although  the  print  outs  use 
blocks).   A  digital  "record"  is  limited  in  the  number  of 
samples  it  can  contain  because  of  Hybrid  computer  memory 
limitations.   This  limit  is  that  the  number  of  samples  t  imes 
the  number  of  channels  simultaneously  digitized  must  be 

equal  to  or  less  than  6144  channel-samples .   As  an  example, 

record 
consider  the  following  hypothetical  problem: 
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Analog  record  length  «*  3,000  sees 

Number  of  analog  channels  containing  data  =  5  channels 

Sampling  rate  =   100  samp les 

sec 

Number  of  channels  digitized 

simultaneously      =  5  channels 

Samples  per  record  selected  =  1024  samples 

record 


1024 — : —  x  5  channels  =  5120  channel-sample    ,,., 

record  ; ' <  6144 


,._.  samples    1  sec 
1024  - — ; —  x 


record     100  samples 


=  10.24 


record 

seconds 
record 


Digital  record  length  =  10.24  seconds 

record 

Note  that  the  value  for  samples  for  record  was  chosen  as 

an  integral  power  of  two  ,  i.e.  2   ,   The  reason  is  that  the 

maximum   number  of  data  points  (samples)  that  can  be  stored 

in  the  IBM  360/67  computer  memory  under  the  FTOR  program  are 

13 
8192  or  2    samples.   Therefore,  for  efficiency  purposes  t  ie 

FTOR  program  was  designed  to  read  in  samples  per  record  wl   .ch 
contain  integral  powers  of  two. 

One  other  commonly  used  word  is  "file."   A  "file"  is  a 
finite  length  of  digital  data  containing  any  number  of 
records  from  one  on  up.   The  length  of  a  "file"  is  predeter- 
mined during  the  digitization  phase  and  is  based  on  the 
number  of  records  desired  per  "file."   An  "end  of  file" 
mark  is  placed  on  the  seven  track  tape  after  the  desired 
number  of  records  have  been  computed. 
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